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VARIATIONAL  methods  of  convolution  integral  and  of 
LARGE  spring  CONSTANTS  -  A  NUMERICAL  COMPARISON 

Julian  J.  Wu 

U.S.  Army  Armament  Research  and  Development  Command 
Benet  Weapons  Laboratory,  LCWSL 
Watervliet,  NY  12189 


nilf.  SEMMAHL.  F^ite  element  solution  formulations  have  been  carried 
Mn  t  !!****  e  value  problem  based  on  two  different  varla- 

®tai“ents:  that  of  convolutional  integral  developed  by  Gurtin 
ana  tnat  of  large  spring  constants  adapted  by  this  writer  for  initial 
N"merlcal  results  indicate  that  both  generate  conver- 
®  J  1 °  the  glven  initial  value  problem  of  a  spring-mass 

system  subjected  to  a  harmonic  forcing  function. 

ii — INTRODUCTION.  Through  a  simple  initial  value  problem,  this 
no  e  demonstrates  the  use  of  the  finite  element  discretization  in  con- 

result  s”  ”12  two  variational  formulations  and  compares  the  numerical 
results.  The  variational  principles  of  convolutional  integral  for 

a^Pr°«emS  were  developed  by  Gurtin  some  sixteen  years  ago  (ref. 
„Since  ^en  these  formulations  have  been  applied  to  obtain 
tranS^ent  probleras  (ref.  3  and  4)  .  However,  the  time  dim- 
treated  separately  from  the  spatial  dimensions  in  the  finite 
*r*“e*5  aPPr°ximation  schemes.  The  viewpoint  adopted  in  this  note  is 
pcca  „  ®fparate  treatment  of  spatial  and  time  coordinates  is  unnec¬ 
essary.  Since  the  initial  value  problems  are  nonself-adjoint,  the 

adlMn?0??1^  var*abional  problems  can  be  formulated  with  the  help  of 

sow? 1  iVaria  l!S  and  thuS  Can  be  USed  ln  PItz-f inite  element 

S**  SUC  formulation  is  used  here  to  compare  with  the  for- 

a  slmnl?  convolution  integral  in  terms  of  numerical  results  for 

a  simple  initial  value  problem. 

equation  nf  slmple  ™afs-spring  system.  The  differential 

equation  of  the  displacement  u(t),  a  function  of  time  t  is 

mu  +  ku  =  fQ  cos  oift  (1) 

where  m  is  the  mass,  k,  the  spring  constant.  A  dot  (•)  denotes  differ- 
entiation  with  respect  to  t.  The  parameters  fQ  and  wf  denote  magnitude 

and  frequency  respectively,  of  the  forcing  function.  The  initial  condi- 
tions  are  given  as 


u(0)  =  Uq  ,  u(0)  =*  u^ 


(2) 
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We  shall  further  use  the  equation 


Thus  Eq,  (1)  has  the  form 

u  +  I02u  “  f  cos  (3) 

2.  VARIATIONAL  FORMULATION  OF  CONVOLUTIONAL  INTEGRALS.  The  varia¬ 
tional  principle  for  the  problem  defined  by  Eqs.  (3)  and  (2)  is  (ref.  2): 

<5l(u)  -  0  (4a) 

where 

I  “  \  [u(t)*u(t)  +  w2t*u(t)*u(t) ] 

"  tu0  +  ~T  +  uj_t  -  cos  uft]*u(t)  (4b) 

The  operator  *  defines  a  convolution  integral  in  the  following  equation 

.t 

u(t)*v(t)  =  J  u ( t— T ) v (t ) dx  (5) 

o 

Vhere  u(t)  and  v(t)  are  two  arbitrary  functions  of  t. 


To  see  that  the  variational  problem  of  Eqs.  (4)  is  indeed  equiva¬ 
lent  to  the  original  problem  defined  by  Eqs.  (3)  and  (2),  one  writes, 
from  Eqs.  (4): 


61  -  [u(t)  +  w2t*u(t)  -  (u0 
for  arbitrary  <5u(t).  Thus,  61 


+  u^t 


cos  toft)]*6u(t) 


“  0  leads  to  Eq.  (6) 


0 


u(t)  +  tl)2t*u(t)  -  (urt  +  -L.  +  n  t  -  JL  COS  U)  t)  -  0  (6) 

u  <i)f  1  U)f  f 

It  is  clear  from  Eq.  (6)  that  u(0)  =  uQ. 

Differentiate  Eq.  (6)  once,  one  has 

u(t)  +  io2  f  u(x)dt  -  u,  -  —  cos  u.t  “  0  (7) 

o  1  »f  f 

Ecl»  (7)  gives  u(0)  =  u£.  Thus  both  of  the  initial  conditions  are  sat¬ 
isfied.  The  differential  equation  is  recovered  when  Eq.  (7)  is  differ¬ 
entiated  once  more.  Note  that  in  obtaining  Eq.  (7)  the  following 
differentiation  formula  has  been  used. 
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Let 


t 

F(t)  “  /  v ( t— T ) u ( T ) dT 
o 

Then 

(t-T)u(T)dT  +  v(0)u(t) 

2i _ VARIATIONAL  FORMULATION  WITH  A  LARGE  "SPRING11  CONSTANT.  Consider 

the  following  variational  problem 

6l(u,v)  -  0  (8a) 

with  1  i 

I(u,v)  =*  -  /  dt  +  J  (w2u-f)v  dt 
0  0 

+  a[u(0)  -  Uq ] v (1)  -  u.jv(0) 

In  Eqs.  (8),  u(t)  is  the  physical  field  variable  and  v(t)  is  the  adjoint 
variable.  This  variational  problem  is  unconstrained  since  the  trial 
functions  of  neither  u(t)  nor  v(t)  are  subject  to  any  end  condition 
requirements.  To  see  that  the  set  of  Eqs.  (8)  is  equivalent  to  the 
original  initial  value  problem,  it  is  only  necessary  to  carry  out  the 
first  variation  and  perform  once  integration-by-part.  Thus,  one  has 

6l(u,v)  ■  0 

1 

“  /  (u+w2u-f)6v  dt 

0 

+  {a[u(0)  -  uQ]  -  u(l) }6v(l)  +  [i(0)  -  Ul]6v(0) 

+  /  (v-Hjj2v)6u  dt 

0 

+  {a  v (1)  +  v(0) }6u(0)  v(l)6u(l)  (9) 

It  is  clear  then  if  one  chooses  v(t)  =  0  and  let  6(t)  be  completely 
arbitrary,  Eqs.  (9)  reduce  to  the  original  initial  value  problem  as  a 
approach  to  infinity. 

A; PROCESS  OF  FINITE  ELEMENT  DISCRETIZATION.  In  case  of  convo¬ 
lutional  integral,  the  variational  equation  used  is  Eq.  (6)  in  Section 
2.  Rewrite  Eq.  (6)  as 

61  =  [u(t)  +  w2t*u(t)  -  F(t) ]*6u(t)  =  0  (10a) 
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where 

F(t)  -  un  +  — j  +  u.  t - =-  cos  U  t 

0  C0f  1  U)f  f 

In  Eq.  (10a),  there  are  three  convolution  integrals  to  be  evaluated 
(a)  u(t)*6u(t);  (b)  w2t*u(t)*5u(t) ;  and  (c)  F(t)*6u(t). 


(a)  For  u(t)*6u(t): 

u(t)*6u(t)  ** 


Let  - 


T 


t 

/  u(t-x)du(x)dx 

o 


“  x/t 


one  has  then 


Consider 


This  integral 


1  _ 

u(t)*6u(t)  »  t  /  u(l-x)6u(x)dx 
0 


1 

I  -  /  u(l-T)6u(x)dT 
0 

is  evaluated  by  finite  element  discretization. 


I 


L 

l 


i=l 


u(l-x)6u(x)dx 


Let 


Hence 


lo  -  0  •  h  - 1  •  h  -  i 

£  -  ■  L  X  -  i  +  1 

X  =  i  K  +  i  "  1] 

dx  =*  L  d£ 

Li 

u(s)  =  u[L  (C+i-1)]  =  u(1)(0 
u(l-x)  =  u[L  (L-5-i+l)] 

=  G[L  (l-£+L-i+l) ]  =  G(L"i+1)(l-0 


(10b) 
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Thus 


1  "  l  l  !  u(L_i+1)(l-5)6u(i)(Od^ 

i=l  L  0 

Use  the  matrix  representations  for  the  shape  function  and  generalized 
coordinates.  One  writes 


u(i)(S)  =  aT(S)  u(i) 


and 


Thus 


Or 


where 

Hence 


U(L“i+l)(i-C)  -  aT(l-^)U(L_i+1^ 

Ia  I  t  <SU(i)  /  a(5)aT(l-5)d5  u(L_i+1) 

i=l  0  ■” 


I-i[  6U(1^  A  U(L-i+1> 
L  i-1  ~ 


A  -  /  a(?)aT(l-5)d5 

0  ~  ~ 


u(t)*6u(t)  =  tl  -  ~  l  U^A  u^L_i+1^ 


i=>l 


(b)  For  to2t*u(t)*$u(t) ; 


(Ha) 


(lib) 


The  evaluation  of  this  double  convolution  integral  is  somewhat  more 
complicated.  First  consider 


t*u(t)  »  /  (t-x)u(x)dx 


Then 


[t*u(t) ]*6u(t) 
t  t-A 

!  {/  (t-A-x)u(x)dx}6u(A)dA 
o  o 
t  t-A 

/  /  (t-A-x)u(x)6u(A)dxdA 

o  o 
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Again  let 


» 


Thus  u(t)  becomes  u(r)t 


T 


t 


6u(X)  to  6u(X)f  etc. 


One  has 


I  -  ts  /  /  (l-X-x)u(x) 6u(X)6xdX 

0  0 


(12) 


It  should  be  pointed  out  that  the  change  of  variables  from  T,X  to  x,X 
(so  that  the  limit  of  integration  is  changed  from  t  to  unit)  is  carried 
out  after  writing  down  explicitly  the  double  convolutional  integral 
and  not  before.  This  is  due  to  the  fact  that  the  definition  of  a 
convolution  integral  requires  that  t  appears  explicitly  in  the  inte¬ 
grals.  To  evaluate  I  of  Eq.  (12)  we  write 

I  -  tsI 


and  work  on  I  instead. 

1  1-X 

I  -  /  /  (l-X-x)u(x)6u(X)dxdX 

0  0 


The  area  of  integration  in  (X,T)  plane  is  the  triangle  bounded  by  lines 
X  ■*  0,  T  ■  0  and  x  "  1— X  (shown  in  shaded  area  in  Figure  1)  .  Using  the 
step  function 


H(l-X-x)  *  1,  X  <  1  -  X 

0,  x  >  1  -  X 

one  can  write 

11 

Is*/  /  H(l-X-x)  (1-X-x)  u  (x)  <5u  (X)  dxdX  (13) 

0  0 


Equation  (13)  will  be  used  for  finite  element  discretization.  We  shall 
divide  the  unit  square  in  (X,x)  plane  into  smaller  squares  of  L  x  L 
(Figure  2).  Let 


5  -  5^  »  LX  -  i  +  1 
n  -  n(i)  -  Lx  -  j  +  1 


(14) 


6 


(0,1) 


T 


Figure  1.  Area  of  integration  for  a  double  integral 
of  convolution:  t*u(t)*6u(t) . 


Figure  2.  Area  of  integration  using  finite  elements. 
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Thus 


6y(X)  -  6y(i)(0 


y(T) 


-(j) 


(n) 


1-X-T-^{L  +  i-  i-  j  +  (i-s-n)} 


and 


L  L  1  1 


1-77  l  l  J  J  H(iJ)(L+i-i-j+i-5-n)  • 

L  1-1  j-1  0  0 


(L+i-i-j+i-5-n)y(j)  (n)6y (Odnd? 


U) 


Or, 


with 


L  L 

I  -  l  l 

i=»l  j=*l  J 

I±1  -  /  /  H(ij)(L+i-i-j+i-5“n)  • 

0  0 

•  (L+i-i-j+i-^-n)y(j)  (n)6y(i)  (Od£dn 


since  1  -  X  -  T  =  0  ,  L  +  1-1J  +  (l-£-r|)  =  0  . 

Or, 

1-5-n-i  +  J-  a+D 

Thus,  three  cases  to  consider  for 
CD  hW) 


(ii) 


H 


(ij) 


1  ,  1  +  j  <  L  +  1 
0  ,  1  +  j  >  L  +  1 


(ill)  -  H(l-^-n)  ,  1  +  j  “  L  +  1 

For  case  (i),  one  has 

i  -  /  J1  (L+2-i-j)y(J)(06y(i)(Odnd5 
J  0  0 


SY(i)T  J1  J1  (L+2-i-j-5-n)a(5)aT(n)dndC  Y(J) 
0  0  ~ 

-  <SY^ 


(15) 


(16) 


(17) 


(18) 
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A(1J)  -  /  /  (L+2-l-j-5-n)a(5)aT(5)dnd5 
0  0  ~ 


1  1 


For  case  (ii). 


hi  - 0 


For  case  (111), 


1  1 


Xij  "  /Q  !  H(1-^“0)(l-C-ri)y(^(ri)6y(1)(5)dnd5 


1 

"  /  /  (i-5-n)y(J)(ri)6y(i)(5)dnd5 

o  o 

^  I  a(£)  J  (l-^-n)aT(n)dnd^  Y^ 
0  ~  0 


Consequently, 


-  6Y(i)T  A 

•*  Of  Of 

-  1  1-5 

a  -  /  a (5)  /  (i-5 -n) aT(n)dnd5 

Oo  ~ 

.s  L  L 

l  l  hi 

L  i-i  j-i  J 

i  +  j<L  +  l  -*  I  -  6Y(i)A(JJMJ> 

ij  **  ~ 

i  +  j  >  L  +  1  -  l±i  -  0 
1  +  j-  L  +  l  1^  -  <SY(1^A  Y(J) 


XUj) 


1  1 


/  /  (Lf2-l-j-5-n)a(5)aT(n)d5dri 
0  0  ~ 


-  1  1-5 

A  -  /  1^0  /  U-5-n)aT(o)dnd5 
0  0 


(19) 


(20) 


(21) 
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And  thus 


2  a  L  L 
art 


u>2t*u(t)*6u(t)  -  w2t*  I  -  l  l  I 

L  i=l  J-l  1J 

(c)  For  F(t)*6u(t)  with  F(t)  given  in  Eq.  (10b),  one  has 
F(t)*6u(t)  ■  (a+bt+c  cos  u^t)*6u(t) 

■  a[l*6u(t)]  +  b[t*6u(t)]  +  c [cos  Wft]*6u(t) 
where,  from  Eq.  (10b): 


(22) 


(23) 


a  "  uo  +  — 5  »  b  "  u. 

0)  ^ 
f 


c  ■  - 


U). 


Now,  for  Eq.  (23),  one  has 

L  T  l 

l*6u(t)  =  7  l  /  a(5)d? 

i=l  0 


2  b  /.xj  1  1 

t*6u(t)  =  ff  l  <$UU'  { (L-i+1)  /  a(Od£  -  /  5a(5)d£} 


i-1 


(24) 


(25) 


(26) 


and 

t  r  (i)T  \  Od«t 

cos  wft  *<$u(t)  =7-1  6UW  J  a(0  cos  [~~  (L-i+l-5)]d5  (27) 

L  i»l  ~  o  ~  L 

Now,  Eqs.  (lib),  (22),  (23)  through  (27),  a  global  matrix  equation  can 
be  written  as 


Or 


<$UT  K  D  ■  6UT  F 

•V  >SI  ^ 

(28) 

KU  =  F 

(29) 

which  is  then  solved. 


The  finite  element  discretization  procedure  for  the  variational  for¬ 
mulation  using  a  large  spring  constant  has  been  described  elsewhere  (see, 
for  example,  ref.  5)  and  will  not  be  repeated  here. 
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5.  NUMERICAL  RESULTS.  Numerical  values  of  the  parameters  in  the 
given  example  as  stated  in  Section  1  are  as  the  following: 

m  ■  1.0  ,  k  ■  1.0  ,  f  ■  1.0  »  u)^  m  0.5 
y0  "  1*°  *  yl  "  1'° 

Computational  results  are  presented  in  Tables  1  through  4.  Table  1  and 
2  compare  results  of  the  two  methods  in  an  interval  of  0  <_  t  <  10,  which 
is  about  the  time  for  a  complete  forcing  cycle.  The  results  for  y(t) 
and  y (t)  are  excellent  for  both  methods.  As  the  interval  becomes  shorter, 
0  jS  t  ^  2  as  shown  in  Table  3  and  4,  the  convergence  is  further  improved. 

TABLE  1.  NUMERICAL  COMPARISONS  BETWEEN  TOO 
UNCONSTRAINED  VARIATIONAL  METHODS 


0_A  t  ^  10.0  10  Elements 


^xXt) 

Convo 

Spring 

Exact 

t  \ 

Integ.  M 

Const.  M 

Solution 

0 

0.999 

1.000 

1.000 

2.0 

1.769 

1.770 

1.768 

4.0 

-1.094 

-1.094 

-1.094 

6.0 

-1.920 

-1.920 

-1.919 

8.0 

0.167 

0.167 

0.167 

10.0 

0.113 

0.114 

0.114 

TABLE  2.  NUMERICAL  COMPARISONS  BETWEEN  TOO 

UNCONSTRAINED 

VARIATIONAL  METHODS 

Oitil.o 

10  Elements 

vy'(t) 

Convo 

Spring 

Exact 

jl/Ss 

Integ.  M 

Const.  M 

Solution 

0 

1.011 

1.004 

1.000 

2.0 

-0.675 

-0.675 

-0.674 

4.0 

-1.520 

-1.518 

-1.512 

6.0 

0.780 

0.778 

0.773 

8.0 

0.691 

0.690 

0.689 

10.0 

-0.391 

-0.385 

-0.381 
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TABLE  3.  NUMERICAL  COMPARISONS  BETWEEN  TWO 
UNCONSTRAINED  VARIATIONAL  METHODS 


0  £  t  £  2.0 


7*F 

Convo 

Spring 

Exact 

t  X 

Integ.  M 

Const,  M 

Solution 

0 

1.000000 

1.000000 

1.000000 

0.4 

1.389154 

1.389154 

1.389153 

0.8 

1.713203 

1.713203 

1.713203 

1.2 

1.911703 

1.911702 

1.911701 

1 

1.938251 

1.938251 

1.938249 

2.0 

1.768413 

1.768416 

1.768416 

TABLE  4.  NUMERICAL 

COMPARISONS  BETWEEN  TWO 

UNCONSTRAINED  VARIATIONAL  METHODS 

0  £  t  i  2.0 

10  Elements 

y^Ct) 

Convo 

Spring 

Exact 

t  \ 

Integ.  M 

Const.  M 

Solution 

0 

0.99999 

1.00000 

1.00000 

0.4 

0.91844 

0.91843 

0.91842 

0.8 

0.67623 

0.67622 

0.67621 

1.2 

0.29662 

0.29662 

0.29661 

1.6 

-0.17425 

-0.71424 

-0.17425 

2.0 

-0.67425 

-0.67413 

-0.67403 

In  conclusion,  we  have  observed  that  the  numerical  convergence  of  the 
method  of  large  spring  constants,  in  the  simple  example  given,  is  at  least 
as  good  as  that  of  the  formulation  through  the  variational  principle  of 
convolutional  integrals.  Both  are  easily  adapted  for  finite  element 
discretization.  Due  to  the  fact  that  the  variational  principles  of  con¬ 
volutional  integrals  can  be  formulated  only  for  a  very  restricted  class  of 
problems  (of  constant  coefficients,  for  example).  The  alternate  approach 
large  spring  constants  appears  to  be  quite  attractive  to  obtain  solu¬ 
tions  of  non-self-adjoint  problems  in  general  and  of  initial  value  prob¬ 
lems  and  initial  boundary  value  problems  in  particular. 


12 


REFERENCES 


1*  H,  E.  Curtin,  "Variational  Principles  for  Linear  Elastodynamics," 
Archive  for  Rational  Mechanics  and  Analysis.  16  (1),  p.  34,  1964. 

2.  H.  E.  Gurtin,  "Variational  Principles  for  Linear  Initial-Value  Prob¬ 
lems,  Quarterly  of  Applied  Mathematics,  22  (3),  p.  352,  1964. 

3.  R.  E.  Nickell  and  J.  L.  Sackraan,  "The  Extended  Ritz  Method  Applied 
to  Transient  Coupled  Thermoelastic  Boundary  Value  Problems,"  SEL 
Report  67-3,  University  of  California,  Berkeley,  CA,  1967. 

4.  J,  Ghahoussi  and  E,  L,  Wilson,  "Variational  Formulation  of  Dynamics 
of  Fluid-Saturated  Porous  Elastic  Solids,"  Journal  of  Engineering 
Mechanics  Division,  Proceedings  of  the  American  Society  of  Civil 
Engineers,  EM4,  p.  947,  1972. 

5.  J.  J.  Wu,  "Solutions  to  Initial  Value  Problems  by  Use  of  Finite 
Element-Unconstrained  Variational  Formulations,"  Journal  of  Sound 
and  Vibration,  Vol.  53,  1977,  p.  341-356. 


13 


TECHNICAL  REPORT  INTERNAL  DISTRIBUTION  LIST 

NO.  OF 
COPIES 

CCMMANDER  1 

CHIEF,  DEVELOPMENT  ENGINEERING  BRANCH  1 

.  ATTN:  DRDAR-LCB-PA  ! 

-DM 

-DP  l 

-DR  ! 

-DS 

-DC  ! 

CHIEF,  ENGINEERING  SUPPORT  BRANCH  1 

ATTN:  DRDAR-LCB-SE  ! 

-SA  2 

CHIEF,  RESEARCH  BRANCH  2 

ATTN:  DRDAR-LCB-RA  1 

-RC  i 

-RM  i 

-RP  1 

CHIEF,  LWC  MORTAR  SYS.  OFC .  1 

ATTN:  DRDAR-  LCB-M  ! 

CHIEF,  IMP.  SIMM  MORTAR  OFC.  1 

ATTN:  DEDAR-LCB-I  x 

TECHNICAL  LIBRARY  •  .  5 

ATTN:  DRDAR- LCB-TL 

TECHNICAL  PUBLICATIONS  &  EDITING  UNIT  2 

ATTN:  DRDAR- LCB-TL 

DIRECTOR,  OPERATIONS  DIRECTORATE  1 

DIRECTOR,  PROCUREMENT  DIRECTORATE  1 

DIRECTOR,  PRODUCE  ASSURANCE  DIRECTORATE  1 


NOTE:  PLEASE  NOTIFY  ASSOC.  DIRECTOR,  BENET  WEAPONS  LABORATORY,  ATTN: 
DRDAR- LCB-TL,  OF  ANY  REQUIRED  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST 


NO.  OF  NO.  OF 

COPIES  COPIES 


ASST  SEC  OF  THE  ARMY 
RESEARCH  &  DEVELOPMENT 
ATTN:  DEP  I  OR  SCI  &  TECH 
THE  -PENTAGON 
WASHINGTON,  D.C.  20315 

COMMANDER 

US  ARMY  MAT  DEV  &  READ.  COMD 
ATTN:  DRCDE 
5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 

COMMANDER 

IIS  ARMY  A RRADC CM 

ATTN:  DRDAR-LC 

-1CA  (PTASTICS  TECH 
EVAL  CEN) 

-LCE 

-ICM 

-ICS 

-ICW 

-TSS(STINFO) 

DOVER,  NJ  07801 

COMMANDER 
US  ARMY  ARRCCM 
ATTN:  DRSAR-LEP-L 
ROCK  ISLAND  ARSENAL 
ROCK  ISLAND,  IL  61299 

DIRECTOR 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRDAR-TSB-S  (STINFO) 

ABERDEEN  PROVING  GROUND,  MD  21005 

COMMANDER 

US  ARMY  ELECTRONICS  CCMD 
ATTN:  TECH  LIB 
FT  MONMOUTH,  NJ  07703 

COMMANDER 

US  ARMY  MOBILITY  EQUIP  R&D  CCMD 
ATTN:  TECH  LIB 
FT  BELVOIR,  VA  220  bO 


COMMANDER 

US  ARMY  TANK-AUTMV  R&D  CCMD 
1  ATTN:  TECH  LIB  -  DRDTA-UL 
MAT  LAB  -  DRDTA-RK 
WARREN  MICHIGAN  4-8090 

COMMANDER 

IJS  MILITARY  ACADFMY 
1  ATTN:  CHMN,  MECH  EN  U'  KPT 
WEST  POINT,  NY  10996 

COMMANDER 
REDSTONE  ARSENAL 
ATTN:  DRSMT-HB 
1  -RRS 

1  -  RSM 

ALABAMA  35509 

1 

1  C  CMMANDER 

1  ROCK  ISLAND  ARSENAL 

1  ATTN:  SARRI-ENM  (MAT  SCI  DIV) 

2  ROCK  ISLAND,  IL  61202 

COMMANDER 

HQ,  US  ARMY'  AVN  SCH 
ATTN:  OFC  OF  THE  LIBRARIAN 
1  FT  RUCKER,  ALABAMA  36362 

C  CMMANDER 

US  ARMY  FGN  SCIENCE  &  TECH  CEN 

ATTN:  DRXST-SD 

220  7TH  STREET,  N.E. 

1  CHARLOTTESVILLE,  VA  22901 

C  CMMANDER 

IJS  ARMY  MATERIALS  &  MECHANICS 
RESEARCH  CENTER 

1  ATTN:  TECH  LIB  -  DRXMR-PL 
WATERTOWN,  MASS  02172 


1 


NOTE:  PLEASE  NOTIFY  C CMMANDER,  A RRADC CM,  ATTN:  BENET  WEAPONS  LABORATORY, 

DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


1 

1 


1 


2 

1 

1 


1 


1 


1 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT) 


NO.  OF 
COPIES 


COMMANDER 

US  ARMY  RESEARCH  OFFICE- 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
COMMANDER 

US  ARMY  HARRY  DIAMOND  LAB 
ATTN:  TECH  LIB 
2,100  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACT 

ATTN:  DRXPE-MT 

ROCK  ISLAND,  IL  61201 

CHIEF,  MATERIALS  BRANCH 
US  ARMY  R&S  GROUP,  EUR 
BOX  65,  FPO  N.Y.  09510 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 
ATTN:  CHIEF,  MAT  SCIENCE  DIV 
DAHLOREN,  VA  22448 

DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  DIR,  MECH  DIV 

CODE  26-27  (DOC  LIB) 
WASHINGTON,  D.  C.  20375 

NASA  SCIENTIFIC  &  TECH  INFO  FAC 
P„  0.  BOX  8757,  ATTN:  ACQ  BR 

raltimore/washingtcn  intl  airport 

MARYLAND  21240 


COMMANDER 

DEFENSE  TECHNICAL  INFO  CENTER 
1  ATTN:  DTIA-TCA 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 

METALS  &  CERAMICS  INFO  CEN 
1  RATTELLE  COLUMBUS  LAB 
505  KING  AVE 
COLUMBUS,  OHIO  4  3  201 

MECHANICAL  PROPERTIES  LATA  CTR 
BATTELLE  COLUMBUS  LAB 
1  505  KING  AVE 

COLUMBUS,  OHIO  43201 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 
1  ATTN:  DRXSY-MP 

ABERDEEN  PROVING  GROUND 
MARYLAND  21005 


1 


1 

1 


1 


NOTE:  PLEASE  NOTIFY  COMMANDER,  ARRADC CM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


NO.  OF 
COPIES 


12 


1 


1 


1 


